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ABSTRACT 



Context. The blazar CTA 102 (z=1.037) underwent a historical radio outburst in April 2006. This event offered a unique chance to 
study the physical properties of the jet. 

Aims. We used multifrequency radio and mm observations to analyze the evolution of the spectral parameters during the flare as a test 
of the shock-in-jet model under these extreme conditions. 

Methods. For the analysis of the flare we took into account that the flaring spectrum is superimposed on a quiescent spectrum. 
We reconstructed the latter from archival data and fitted a synchrotron self-absorbed distribution of emission. The uncertainties of 
the derived spectral parameters were calculated using Monte Carlo simulations. The spectral evolution is modeled by the shock-in-jet 
model, and the derived results are discussed in the context of a geometrical model (varying viewing angle) and shock-shock interaction 
Results. The evolution of the flare in the turnover frequency-turnover flux density (v„,-S„,) plane shows a double peak structure. The 
nature of this evolution is dicussed in the frame of shock-in-jet models. We discard the generation of the double peak structure in the 
ym-Sm plane purely based on geometrical changes (variation of the Doppler factor). The detailed modeling of the spectral evolution 
favors a shock-shock interaction as a possible physical mechanism behind the deviations from the standard shock-in-jet model. 



Key words, galaxies: active, 
individual: CTA 102 



1. Introduction 



- galaxies: jets, - radio continuum: galaxies, - radiation mechanisms: non-thermal, - galaxies: quasars: 



The blazar CTA 102 (B22 30-H114) has a redshift of z = 1.037 
(iHewitt & Burbidgd Il989t) . It is classified as a highly polar- 
ized quasar (HPO) with a linear optical polarization above 3% 
dvir on-Cettv & Veron 2003) and a V-band mag nitud e of 17.33. 
The source was ob s erve d for the first time by Harris & Roberts! 
d 1 960l) . iKardashevI d 1 964) reported on possible si gnals from an 
extrate rrestrial civiUzation coming from CTA 102. Sholomitskiil 
(Il965h found the first variation in flux density for a radio source. 
Later observations identified CTA 102 as a quasar. 

Since that time CTA 102 has been the target of numerous 
observations at different wavelengths. Besides the aforemen- 
tioned variation in the radio flux density, CTA 102 also shows 
variation in the optical band. iPica et al" ] (Il988h reported a vari- 
ation range of 1.14mag around an average value of 17.66 mag 
in 14 years, and an increase of 1.04 mag within two days in 
1978, which is so far the most significant outburst. CTA 102 has 
been monitored since 1986 within the cm-observations of the 
Metsahovi telescope. The strongest radio flare since the begin- 
ning of the monitoring took place around 1997, and a nearly si- 
multaneous outburst in the optical R-band w as observed with 
the Nordic Optical Telescope on La Palma (iTornikoski et al.l 
Ii999). The source has been detected in the y-ray regime by 
the telescopes CGRO/EGRET and FermijEP^ with a luminos- 
ity, Ly - 5 y. lO'^^erg/s, defining CTA 102 as a y-bright source 



(iNolan et alJll"99llAbdo et"ani2009l) . Meriin and VLA observa- 
tions at 2 GHz, 5 GHz and 15 GHz revealed the kpc-scale struc- 
ture of CTA 102 which con sists of a central core and two faint 
lobes. dSpencer et al.lll989h . The brighter lobe has a flux density 
of 170 mJy at a distance of 1.6 arcsec from the core at position 
angle (P. A.) of 143° (measured from north through east). The 
other lobe, with a flux density of 75 mJy, is located 1 arcsec from 
the center at P. A. -43°. The spectral indices between 2 GHz and 
5 GHz of the lobes are -0.7 for the bright and -0.3 for the other 
one. 

High-resolution VLBI observations at 1.4 GHz and 5 GHz 
resolved the central object into three components and a diffuse 
tail bending to the southeast. These observations provide and 
upper limit around 10 c (0.5mas/yr) f or the superluminal mo - 
tion of the components (Baath 1988: IWehrle & Cohenl[T989h . 
Several observations at different frequencies (for example at 
326 MHz) confirmed the el ongation of the source to the south- 
east ( Altschuler et al.'l995'). 

Within t he framework of the VLBA 2 cm-Survey (e.g., 
IZensus et al.l l2002) and the MOJAV^] program (Monitoring of 
Jets in A ctive galactic nuclei with VLBA Experiments) (e.g., 
iLister et al. 2009a) CTA 102 has been monitored regularly, be- 
ginning in mid 1995. These intensive observations deliver a de- 
tailed picture on the morphology and kinematics of this source at 



http://www.physics.purdue.edu/MOJAVE 
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Fig. 1. Radio-mm light curves for CTA 102, centered around the 2006 radio flare. 



15 GHz. Kinematic analysis show apparen t velocities of the f ea- 
tures in the jet between 0.7 c and 15.40 c (iListeret al.ll2009bl) . A 
multifrequency VLBI study i ncluding data at 9 GHz , 43 GHz 
and 22 GHz was reported by 'Rantakvro et al.l (l2003l) . The re- 
sults from the multifrequency VLBI observations were com- 
bined with the continuum monitoring performed at single-dish 
observatories at 22 GHz, 37 GHz, 90 GHz and 230 GHz. Within 
this multifrequency data set (November 1992 until June 1998), a 
major flare in CTA 102 around 1997 was confirmed. The authors 
could conclude that this event was connected to the eject i on of a 
new jet featu r e. The sa me was noted by Savolai nen et al.l(l2002h . 
iJorstad et al"](l2005h and lHovatta et al.l (120091) found Lorentz fac- 
tors, r, of 17 and 15, respectively, and Doppler factors, 6, be- 
tween 15 and 22 associated to this ejection. The 2006 radio flare 
in CTA 102 has been observed at cm - mm total flu x density and 
multifrequency VLBI observations (lFrommll2009h . The results 
of the the multifrequency VLBI observations will be presented 
in second paper 

In this work, we concentrate on the analysis of the cm - mm 
light curves. The organization of the paper is the following. 
In Sect. |2] we present the radio/mm light curves during the 
2006 flare and perform the spectral analysis. The theoretical 
background and the fitting technique are introduced in Sect. [3] 
The results of this analysis are shown in Sect. |4] and are applied 
to the 2006 flare in CTA 102 in Sect.|5] In Sect.|6]we discuss the 
different models that can explain the observations. 

Throughout this paper we define the optically thick spectral 
index as a, > and the convention used for the optically thin 
spectral index is 5 oc v^", where ao < 0. The latter can be de- 
rived from the spectral slope, s, via the relation oq - -{s - l)/2. 



Table 1. Average time sampling and average flux density uncer- 
tainties for the used light curves 



V [GHz] 


Observatory 


(fobs) [days] 


<A5ote> [Jy] 


4.8 


UMRAO 


46 


0.09 


8.0 


UMRAO 


38 


0.07 


14.5 


UMRAO 


32 


0.07 


37 


Metsahovi 


16 


0.24 


230 


SMA 


27 


0.18 


340 


SMA 


132 


0.25 



2. Observations: cm-mm light curves 

For our analysis we focused on the radio flare around April 
2006 and used observations spanning from 4.8 GHz to 340 GHz 
(see Fig. [TJ. The observations have been carried out by the 
Radio Observatory of the University of Michigan (UMRAO), 
the Metsahovi Radio Observatory, and the Submillimeter Array 
(SMA). The average sampling time intervals and flux density 
uncertainties are presented in Table [1] 

Figure [1] shows the total flux densities measured by the tele- 
scopes at different frequencies. The most prominent feature in 
the light curve is the major flare around 2006.2, best seen at 
37 GHz. This feature is surrounded by smaller flares in 2005.2, 
2007.6, 2008.5, and 2009.4. The light curves show the typi- 
cal evolution of a flare: the flaring phenomenon usually starts 
at high frequencies and propagates to lower frequencies with 
a ceitain time delay of the peak, but there are also flares that 
develop simultaneously over a wide frequency range. The flare 
around 2005.0 appears nearly simultaneously at the highest fre- 
quencies (230 GHz and 340 GHz) and delayed at 37 GHz and 
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14.5 GHz, whereas it seems that the lowest frequencies (4.8 GHz 
and 8 GHz) are not affected by this event. 

The main flare shows the typical evolution: it is clearly visi- 
ble at all frequencies with increasing time delays towards lower 
frequencies. Due to the poor time sampling of the 340 GHz ob- 
servations, the rise and the time shift at this frequency are not 
easy to determine. The remarkable double peak structure of the 
230 GHz measurements is an interesting feature and will be dis- 
cussed later After the flare, the flux decreases at all frequencies, 
with a steeper descent at higher frequencies. 

2.1. Time sampling and interpolation 

To perform a spectral analysis of the light curves, simultaneous 
data points are needed. This was achieved by performing a linear 
interpolation between the flux density values from the observa- 
tions. The choice of an adequate time sampling, Af, depends on 
the time interval of the observations and the significance of the 
frequency for the determination of the peak frequency, i.e., the 
frequency where the spectral shape changes from optically thick 
to optically thin. If a too short time interval is chosen, most of the 
data points would be interpolated ones. On the opposite situa- 
tion, if a too long time interval is chosen, most of the light curves 
would be smoothed. Moreover, in this case, the data are sampled 
inhomogeneously at the different frequency, so a compromise 
is required. The influence of a certain frequency on the calcu- 
lation of v,„ can be estimated from the light curve (see Fig. [TJ. 
The turnover frequency is usually between the frequencies of 
the highest and second highest flux density at a certain epoch. 
In the case studied here, the dominating frequencies are 37 GHz 
and 230 GHz. From this distribution it could be concluded that 



an appropriate time sampling should not be significantly shorter 
than the observational time interval used for these frequencies. 
A time sampling of Af = 0.05 yr was selected for the interpola- 
tion, corresponding approximately to the observation cadence of 
the 37 GHz-light curve. 

The correct handling of the uncertainties in the interpolated 
flux densities requires knowledge of a mathematical relation 
describing the light curves. Since this approach is out of the 
scope of this paper, we assigned the maximum uncertainty of 
the two closest observed flux densities to the interpolated flux 
density. For our analysis we focused on the time interval be- 
tween 2005.6 and 2006.8. The interpolation of the poorly sam- 
pled 340 GHz light curve after 2006.2 could induce artificially 
flat spectra. We took this fact into consideration by excluding 
the interpolated 340 GHz flux densities from the spectral analy- 
sis for 2006.2 < t < 2006.8. 

2.2. Spectral analysis 

Following (Tiirler et al.l2000l) . a synchrotron self absorbed spec- 
trum is described by 



S V — S nil 



l-exp(-r„,(y/v„,r-"') 
1 - exp (-T„,) 



(1) 



where t,„ ^ 3/2 ^ -^1 - ^f^ - 1 j is the optical depth at the 

turnover frequency, S „ is the turnover flux density, v,„ is the 
turnover frequency and a, and q-q are the spectral indices for the 
optically thick and optically thin parts of the spectrum, respec- 
tively. 
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The observed spectra may in general be thought as the su- 
perposition of the emission from the steady state and perturbed 
(shocked) jet. We re-constructed this quiescent spectrum from 
archival data and we identified the quiet state with the mini- 
mum flux density of the low frequency light curves (4.8 GHz 
- 37 GHz) around f = 1989.0 (see Fig.|5]and Tabled. The flux 
densities were fitted by a power law S (v) 
tained c = 7.43 + 0.65 and ao = -0.45 + 0.04. 



CgV"° and we ob- 



Table 2. Frequencies and flux density values for the quiescent 
spectrum 



V [GHz] 


S [Jy] 


Observatory 


4.8 


3.68±0.05 


UMRAO 


8.0 


2.95±0.06 


UMRAO 


14.5 


2.10±0.06 


UMRAO 


22 


2.12±0.14 


Metsahovi 


37 


1.62±0.11 


Metsahovi 



For the spectral analysis, we removed the contribution of the 
quiescent spectrum from the interpolated data points and applied 
Eq. [1] The uncertainties of the remaining flaring spectrum were 
calculated using the eiTors of the interpolated data points and 
the obtained uncertainties of the quiescent spectrum. During the 
fitting process we allowed both spectral indices, a, and ao to 
vary. Figure |3] shows the result of a spectral fitting at a selected 
epoch applied to interpolated light curve data. 



2.3. Error analysis of the spectra 

Due to the non-linear nature of Eq. [T]we applied a Monte-Carlo 
simulation to the observed/interpolated flux densities in order to 
derive estimates for the uncertainties of the fitting parameters 
(Sm, Vm, <^t, and ao)- Random values for each simulated spec- 
trum were drawn from a Gaussian distribution with a mean, //, 
equal to the observed flux density, S (v), and a variance, cr, equal 
to the observed flux density eiTor, A5(v). In order to achieve a 
good statistical ensemble, up to 1000 spectra were simulated per 
interpolated spectrum. Each of these simulated spectra were fit- 
ted with a synchrotron spectrum (see Eq. [1}. As the expected 
value and uncertainty of a spectral parameter we take the mean 
and standard deviation of its simulated probability density dis- 
tribution, respectively. The derived spectral parameters and their 
uncertainties are presented in Table |3] 



3. Synchrotron radiation and the shock-in-jet model 

In this section we review the shock-in-jet model as derived by 
lMarscher&Gea3 (fl985) and different modifications to it. In this 
review we have included the redshift dependence in the different 
equations where it is required. In this way, we can fit the model 
to observational data from any source. 



3.1. Synchrotron radiation 

The specific intensity from a power law distribution of relativis- 
tic electrons, N{E) - KE \ where K is the normalization coef- 
ficient of the distribution and s the spectral slope is defined by: 



10^ - 



total 
flare 
quiet 
data 




Flaring Spectrum 

S„, ;7.67 |Jy] 1/ 
a a — 



^(1 



(2) 



where Cy and Ky are the emission and absorption coeflicients, re- 
spectively, and Ty is the optical depth. The emission and absorp- 
tion c oefficients can be written as (for details see iPacholczvfl 
II970h : 



Ky oc KB 



,(i+2)/2^-(s+4)/2 



(3) 
(4) 



where B is the magnetic field and v the frequency in the source 
frame. With this definition of the emission and absorption coef- 
ficients the specific intensity is: 



(1 



(5) 



Depending on the optical depth, Ty, Eq.|5] describes an optically 
thick (Ty > 1) or optically thin (Ty < 1) spectrum with their 



characteristic shapes ly oc v^/~ and ly 



^-(i-l)/2 



respectively. 



The optically thin flux density from a jet located at luminos- 
ity distance D^, with radius R, conical geometry (half opening 
angle (f>), constant velocity /3 and viewing angle ■& is given by 



Fig. 3. Result of the spectral fitting to the 2006.20 data. The 
dashed-dotted line coiTesponds to the quiescent spectrum, the 
dashed line to the flaring spectrum, and the solid black line to 
the total spectrum. The values presented indicate the spectral 
turnover of the flaring spectrum. 



(6) 



where Q - nR IDj{\ + z) is the solid angle, and z was the 
redshift. Taking into account the influence of the Doppler fac- 
tor, D = (1 -yScos)?)"' (see e.g.. Lind & Blandfo rd 1985]), 

where T = -y6^) ^ the Lorentz factor, Eq.|6]leads to: 



Sy^{\+Z) 



-(.5-3)/2 j^l^f. ^(.5+I)/2 ^(.5+3)/2^-(i-l)/2_ 



(7) 
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Table 3. Spectral parameters (expected values) and corresponding uncertainties from the Monte-Carlo simulation 



t[yr] 


v,„ [GHz] 


S,„ [Jy] 


a, 




2005.60 


222±99 


3.3+0.3 


0.33+0.14 


-0.27+0.50 


2005.65 


163±27 


5.9+0.5 


0.56+0.11 


-1.80+0.57 


2005.70 


135±25 


5.7+0.8 


0.59+0.15 


-0.90+0.50 


2005.75 


119±19 


6.5+1.1 


0.70+-0.16 


-0.78+0.42 


2005.80 


109±12 


8.4+1.0 


0.78+0.14 


-1.02+0.28 


2005.85 


96±12 


7.7+1.0 


0.82+0.16 


-0.72+0.25 


2005.90 


88±14 


6.7+0.9 


0.81+0.16 


-0.42+0.28 


2005.95 


73±14 


6.1+0.5 


1.0+0.3 


-0.22+0.15 


2006.00 


62±12 


6.2+0.5 


1.2+0.3 


-0.20+0.11 


2006.05 


56+11 


6.5+0.6 


1.3+0.4 


-0.19+0.11 


2006.10 


63+12 


6.5+0.6 


1.1+0.3 


-0.20+0.13 


2006.15 


60±9 


7.2+0.6 


1.0+0.3 


-0.28+0.13 


2006.20 


57±7 


7.9+0.6 


1.1+0.2 


-0.33+0.12 


2006.25 


56±7 


8.5+0.8 


0.93+0.22 


-0.47+0.20 


2006.30 


56±7 


8.4+0.7 


0.94+0.21 


-0.43+0.19 


2006.35 


53±7 


8.4+0.7 


0.97+0.24 


-0.45+0.20 


2006.40 


41±5 


7.2+0.6 


1.0+0.3 


-0.49+0.13 


2006.45 


27±6 


4.5+0.3 


1.1+0.5 


-0.35+0.13 


2006.50 


25±7 


3.8+0.4 


0.89+0.56 


-0.40+0.18 


2006.55 


25±7 


3.4+0.3 


0.37+0.38 


-0.66+0.32 


2006.60 


25±10 


2.9+0.4 


0.28+0.35 


-0.61+0.35 


2006.65 


25±9 


2.6+0.4 


0.24+0.19 


-0.47+0.34 


2006.70 


21±11 


2.3+0.4 


0.26+0.13 


-0.59+0.34 


2006.75 


17±7 


2.0+0.3 


0.22+0.18 


-0.63+0.33 


2006.80 


16±8 


1.8+0.3 


0.24+0.11 


-0.75+0.43 



The turnover frequency, defined as the frequency where the op- 
tical depth is equal to unity, can be expressed as: 

V,„ cc (1 + z)-' (.^/^B(-'+2)/2^(.,+2)/2^2/<-4) 

The turnover flux density S„, (v,„) is then obtained by inserting 
Eq.[8]intoEq.|7J 

5„, cc (1 + z)DlR' (/^^B2.,.3^5^3...7y/(-4) 

The equations above were successfully used to model the steady- 
state emission from quasars, where a decrease in the particle 
density and the magnetic field B along the jet was assumed 
(e.g., |K5nigl 1981). However, the observed flux densities dur- 
ing flares showed a more complicated behavior and could not be 
described by steady-state models. 

3.2. Shock-in-jet model 

The shock-in-jet model of iMarscher & Gearl (1 19851) studied the 
evolution of a traveling shock wave in a steady state jet. During 
the passage of the shock through a steady jet, the relativistic 
particles are swept up at the shock front and gain energy while 
crossing it. In this model, the flaring flux density is assumed to 
be produced by the accelerated particles within a small layer 
of width X behind the shock front. The width of this layer is 
assumed to depend on the dominant cooling process and can be 
approximated by x oc t^ooi, where fcooi is the typical cooling time 
i.e., synchrotron cooling time. 



3.2.1. Evolutionary stages 

Compton losses: Assuming the photon energy density, Mph, is 
higher than the magnetic energy density, Mb - B^/iSn), the in- 
verse Compton scattering is the dominant energy loss mecha- 
nism during the first stage of the flare. The width of the layer 



behind the shock front during this so-called Compton stage, xi , 
is computed to be 

xiccB1/V'/2d'%i. (10) 

An approximation of the photon energy density, Mph, can be ob- 
tained by integrating the emission coefficient, ey, over the opti- 
cally thin regime (v„, < v < v„iax fl which leads to 

Mphcc7^(B3-7^-5)'/^ (11) 

The final expression for the width of the layer behind the shock 
front during the Compton stage can be written as 

XI cc7;-(-'+5)/8^-l5-3(.+ I)/8^I/2^-l/2_ (12) 

Synchrotron losses: Synchrotron losses become more impor- 
tant at the point where the photon energy density, Mph, is compa- 
rable to the magnetic energy density, Mb. The width of the layer 
X2 in the synchrotron stage can be computed by replacing Mph in 
Eq.[Tn]with Mb - l{2>n) and is given by 

X2 oc B-"^D"^v-"\ (13) 

Adiabatic losses: Radiative losses become less important in the 
last stage of the shock evolution. During this final stage the 
losses are dominated by the expansion of the source and the 
width of the layer, xj,, is assumed to be similar to the radius of 
the jet 

X30CR. (14) 

The evolution of the turnover frequencies v„,,, and turnover flux 
densities 5,,,,,, where / indicates the different energy loss stages 

^ IMarscher & Gead ( 119851) included only first order Compton scat- 
tering 
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(l=Compton, 2=synchrotron and 3=adiabatic loss stage), can be 
derived by inserting the expressions for jc, into Eq.[8]and Eq.|9l 



-(j+4)/(j+5) ^-1/4^1/4 £,(.s+3)/(.s+5) 



S„A oc (l+^p+15)/(2.+ 10) ^-2^11/8^1/.^, 



,(3.s+10)/(.s+5) 



> m,l ^ Vi i-'L 

y,n,2 <^ (1 + zr^'^*^'^'^'^[K^B'-'D'^'f"-'^'^ 

S,„,2 - (l+2)(2-15)/(2.+ 10)^^2^2^^552.-5^3.,+10jl/(^+5) 

V„,3 cc (1+z)-'[/?/^B(-2)/2d("2)/2J2/(-4) 



(15) 
(16) 

(17) 
(18) 

(19) 
(20) 



Furthermore iMarscherl (Il990b : Lobanov & ZensusI (1 1999b as- 
sumed that the evolution of /T, B and D could be written as a 
power-law with the jet radius R (notice that in a conical jet the 
distance along the jet L is linearly proportional to the jet radius 
R,L oc R), 



KccR'' BocR-'' Doc R'' . 



(21) 



Then, a relation between the turnover flux density and the 
turnover frequency can be found: 



Vmj (I+Z)"'/?"' 



S^j oc Df(l+zrR\ 
which leads to 



(22) 
(23) 

(24) 



where e, = /l/n/- The exponents fi, rij, pj and qj include the 
dependence with the physical quantities B, N{E), K, and D and 
can be written as: 

„i = -(t+l)/4-d(s + 3)/(s + 5) (25) 

n2 = -[2k + b(s-l) + d(s + 3)]/(s + 5) (26) 

ns = -[2(k-l) + (b + d)(s + 2)]/(s + 4) (27) 

/i = (n-b)/8-d{3s+lO)/{s + 5) (28) 

/2 = 2-[5k + b(2s -5) + d(3s + lO)]/(s + 5) (29) 

/3 = [2s+l3-5k-b(2s + 3)-d{3s + 7)]/(s + 4) (30) 

pi = -(.v + 4)/(^ + 5) (31) 

P2 = -is + 4)/{s + 5) (32) 

P3 = -1 (33) 

qi = (2i-H 15)/(2i+ 10) (34) 

q2 = (2i + 15)/(2i-H 10) (35) 

^3 = 1. (36) 

The typical evolution of a flare in the turnover frequency - 
turnover flux density (v,„ - 5,„) plane can be obtained by in- 
specting the /^-dependence of the turnover frequency, v,„, and 
the turnover flux density, S,„, for a typical set of parameters: 
b=l, s = 2.5, k - 3 (assuming an adiabatic flow, k = 2{s + 2)13) 
and d - Q (assuming constant velocity). These parameters lead 
to a set of exponents «, and //, which are summarized in Table |4] 
During the first stage, where Compton losses are domi- 
nant, the turnover frequency, v,„, decreases with radius, R, while 
the turnover flux density, 5,,,, increases. In the second stage, 
where synchrotron losses are the dominating energy loss mech- 
anism, the turnover frequency continues to decrease while the 
turnover flux density is constant. Both the turnover frequency 
and turnover flux density decrease in the final, adiabatic stage. 



Table 4. Values of the indices used in the generalized spectrum 
expression of Eq. |24l n, and f, for b - \, s - 2.5, k - 3, and 
d = 



stage 




f. 


e 


Compton 


-0.5 


1.3 


-2.6 


Synchrotron 


-1 








Adiabatic 


-1.3 


-0.8 


0.6 



3.2.2. Generalization of tine sliock-in-jet model 

iTiirler et"al1 (l200d) expanded the shock-in-jet model to non- 
conical jets, e.g., the distance along the jet, L, is no longer a 
linearly proportional to the jet radius, R. This modification is 
expressed hy R oc U with (-1 < r < 1), which represents a 
collimating (r < 0) or expanding (r > 0) jet. An important im- 
plication of allowing non-conical jets (r + 1) in the shock-in-jet 
model is the r-dependence of the evolution of the spectral pa- 
rameters (b, k and d). Therefore the evolution of B, K and D 
along the jet is given by: 



K oc L-'", B oc L-''', 



DocL 



-rd 



(37) 



Be sides the modification of the jet geometry, iTiirler et"al] 
(200(f) parametrized the temporal evolution of the turnover fre- 
quency, v„„ and the turnover flux density, 5„„ using the equations 
of superluminal motiorQ 



t = where Vob.s = -; , . „ . 

Vobs (l+z)sm§ 



(38) 



Using the definition of the apparent speed, ySapp - (j6 sin )?)/(! - 
/3 cos -&), the bulk Lorentz factor, F, and the Doppler factor, D, 
the equation above results in 



t 



PcDT 



(39) 



Assuming that yS ~ 1 and & ~ 1 /F «: 1, it can be shown that F oc 
D (Taylor expanding cos 1 /F in the definition of D). Including 
the non-conical jet geometry (R oc U) the equation above can be 
written as 



t oc {\ +z)D-^R^ I '' oc{\+z)RP , 



(40) 



where p - (2dr + 1)1 r. By replacing R in Eqs. l22l and |23] bv 
Eq. |40] the temporal evolution of the turnover frequency and 
turnover flux density is: 



V,„ cc {\+zfPiP-".)IP .f>IP 



S,„ oc D-^\l +z) 



\(qiP-fi)/p . ifi/p 



(41) 
(42) 



The equations above can be further simplified by replacing Di 
(1 + z), which leads to: 



5,„oc(l +^)[(?.-2)p-/;]/p. j//p_ 



3.2.3. Modification of the compton stage 



(43) 



The determination of the distance that the relativistic par- 
ticles travel behind the shock front before losing most of 
their energy, was a crucial parameter in shock-in-jet models. 
Biomsson & Aslaksen (2000) presented a more accurate deriva- 
tion of this distance than previous works, including the possi- 
bility of multiple Compton scattering (in the Thompson regime) 



The time used here was that in the observers frame. 
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during the first rising phase of the flare. In the case of first order 
Compton scattering, the evolution of the first stage was predicted 
to be less abrupt when compared to Marscher& Gear (1985). 
Therefore the existence of a synchrotron stage was no longer 
needed. The evolution of the turnover-flux density, v,„, and the 
turnover frequency, S with distance (or time) is given, within 
this model, by: 



V,„ oc ^-[4(.5+2)+3i(s+4)]/[3(.5+12)] 
S oc ^-[4(^-13)+6fe(i+2)]/[3(i+I2)] 



(44) 
(45) 



assuming no changes in the Doppler factor and constant velocity 
of the shocked particles. Using the same set of parameters as in 
Table|4]leads to rii = -0.9 and /i = 0.3, resulting in a less steep 
slope e\ - /i/wi = -0.3 in the Vm - S'm-plane which is flatter 
than in lMarscher & Gead(ll985h case. 



4. Results 

The spectral evolution of the 2006 radio flare in CTA 102 is pre- 
sented in the turnover frequency - turnover flux density (v,„ -5,„) 
plane, shown in Fig.|4l Figures|5]and|6]display the time evolution 
of the spectral parameters (v,„, 5,,,, ao, and a,). 



10^ 



Co 



10" 



2006.3 2005.1 




























A' 












2006.8 1 





2005.6 



10^ 



10^ 



[GHzl 



Fig. 4. The 2006 radio flare in the turnover frequency - turnover 
flux density plane. The time labels indicate the time evolution 
and the temporal position of local and global extrema in the spec- 
tral evolution. 



In the next paragraphs, we summarize the evolution of 
the event in the turnover frequency — turnover flux den- 
sity plane and compare it with the standard shock-in-jet 
model jM arscher & Gear 1985) and later modifications to it 
dBiornsson & Aslaksenll2000i) . This information is also summa- 
rized in Table m 

The flare started around 2005.6 with a high turnover fre- 
quency (v„, ~ 200 GHz) and a low turnover flux density (S,,, ~ 



3 Jy). During the first 0.2 yr, the turnover flux density, 5,„, in- 
creased, reaching 5,„ ~ 8.4 Jy, while the turnover frequency 
decreased to v„, ~ 110 GHz. The slope of the optically thick 
part of the spectrum, represented by a,, steepened (0.33 to 0.78), 
while the optically thin spectral index, ao, steepened during the 
first 0.05 yr and flattened (-1.8 to -1.02) afterwards. The large 
uncertainties of the optically thin spectral index, ao, are due to 
the lack of data points beyond 340 GHz. The exponential re- 
lation between the turnover flux density, Sm and the turnover 
frequency, v™, {S„, oc v%, see Sect. I3.2l i. led to a value of 
e = -1.21 +0.22. 

As mentioned in Sect. [3] the dominant loss mechanism dur- 
ing the first stage of flare (increasing turnover flux density 
while the turnover frequen c y was decreasing) was Compton 
scattering lMarscher & GeaJ (Il985h p r edict ed a value of e = 
-5/2, whereas iBiornsson & AslaksenI (12000 ) derived e - -0.43 
using a modified expression for the shock width (both assumed 
s = 2.4 and b - I). The obtained value of e = -1.21 + 0.22 was 
in between these two values, but it is impossible t o reproduce 
using the approach of IBiornsson & AslaksenI (|2000|) . 

During the time interval between 2005.8 and 2006.0, in- 
volving 0.2 yr, the turnover flux density and turnover frequency 
decreased to S„, ~ 6.2 Jy and to v„, ~ 62 GHz, respectively. 
The slope in the v,„ - 5,,,-plane changed to e = 0.77 + 0.11. 
The average optically thick spectral index reached a value a, ~ 
0.82 + 0.17 while the optically thin spectral index, ao, continued 
rising to -0.22. This behavior of the turnover values fits well 
in the adiabatic s t age in the shock-in-jet model. For this stage 
iMarscher & Gearl (Il985b derived an exponent e = 0.69 (assum- 
ing s - 3), which is well within our value. 

However, the increase of the turnover flux density starting in 
2006.0, which reached a peak value of S,,, ~ 8.5 Jy in 2006.3, 
cannot be explained within the frame of the shock-in-jet model. 
Its behavior, though, resembles that expected from a Compton 
stage (e = -0.99 + 0.46). The spectral indices during this stage 
were roughly constant, a, = 1.33 + 0.27 and ao = -0.26 ± 0.13. 

After 2006.3, the turnover flux density and turnover fre- 
quency decreased. This last phase could be understood as an adi- 
abatic stage and the power law fit in the v,„ - 5,,,-plane led to an 
e - 1.24 + 0.10. The optically thin spectral index was nearly 
constant, ao = -0.46 + 0.21 while the evolution of the opti- 
cally thick spectral index, a,, could be divided into two parts: 
a, = 0.96 ± 0.28 until 2006.5 and a, = 0.26 + 0.18 after 

Our results show no evidence for a synchrotron stage, char- 
acterized by a nearly constant turnover flux density, Sm, while 
the turnov er frequency, v,„, is decreasing (e = -0.05 for s - 2.4 
and b = 1 . IMarscher & Gear 1985). 

The first hump in the evolution of the flare in the v,„ - Sm- 
plane (2005.6 < f < 2006.0) fulfills, despite no evidence for 
a synchrotron stage, th e phenomenological requ irements of the 
shock-in-jet model of IMarscher & Gead (Il985l) . However, the 
evolution after 2006.0 did not follow the predicted evolution. 
The turnover flux density and the turnover frequency should con- 
tinue decreasing. But their behavior mimicked the evolution of 
a Compton stage (2006.0 < f < 2006.3) and the one of an adi- 
abatic stage (2006.3 < t < 2006.8). Considering these points, 
we decided to model the evolution of the 2006 flare in C TA 102 
according to the standard shock-in-jet model (Marscher & Gearl 
119851) including a second Compton and adiabatic stage. This de- 
cision is evaluated in Sect. |6] 
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Table 5. Different stages of the spectral evolution and their char- 
acteristics 



time 


stage 


model name 


f (5,,, « v^„) 






2005.6 - 2005.8 


Compton 


CI 


-1.21+0.22 


rising 0.33 - 0.78 


rising -1.80- -1.02 


2005.8 - 2006.0 


Adiabatic 


Al 


0.77±0.11 


0.82+0.17 


rising -1.02- -0.21 


2006.0 - 2006.3 


Compton-like 


C2 


-0.99+0.46 


1.03+0.27 


-0.26+0.13 


2006.3 - 2006.8 


Adiabatic 


A2 


1.24+0.10 


0.96+0.29", 0.26+0.18* 


-0.46+0.21 



until 2006.4, * after 2006.45 




2005.6 2005.8 2006.0 2006.2 2006.4 2006.6 2006.8 



: Ml 

1,1 


, 











2005.6 2005.8 2006.0 2006.2 2006.4 2006.6 2006.8 

t [yr] 



Fig. 5. Temporal evolution of the 2006 radio flare: top turnover 
frequency and bottom turnover flux density. The dashed vertical 
lines correspond to the time labels in Figure |4] and indicate the 
extrema in the evolution. 
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2005.6 2005.8 2006.0 2006.2 2006.4 2006.6 2006.8 

f [yr] 
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1.5- 




2005.6 2005.8 2006.0 2006.2 2006.4 2006.6 2006.8 

' [yr] 

Fig. 6. The temporal evolution of the optically thin {ao, top) 
and thick (ff,, bottom) spectral indices for 2006 radio flare. The 
dashed vertical lines corresponds to the time labels in Figure |4] 
and indicate the extrema in the evolution. 



5. Modeling the 2006 radio flare 

In this section we present the results of applying the shock-in-jet 
model and the fitting technique to the observed spectral evolu- 
tion of flie 2006 radio flare in CTA 102. 

5.1. Fitting technique 

We used a multi-dimensional -optimization for deriving a 
set of parameters that fit the spectral evolution of the different 
stages. Our approach consists on fitting the temporal evolution 
of the spectral turnover values, v,„ and S,„ using Eqs. [41] and 
|43]and the definition of the spectral exponents (see Eqs. l25ll36] ). 
The proportionalities in Eqs. 141] and |43] can be removed by in- 
troducing the constants Cy,^ , and cs_ , which reflect logarithmic 
shifts of the turnover frequency and flux density and depend on 
the intrinsic properties of the source and flare, having no further 
importance for our study. From the observed values, v^^**'^ and 

S (i indicating the radiation loss stage, and j indicating the 
position among the total number of points in the stage, q, so that 
j - 1 ■ ■ - q), the proportionality constants can be derived as: 



cs = S -rl''", (47) 

where tobsj is the time in the observers frame. Using these def- 
initions, the constant for the spectral evolution in the v,„ - 5„, 
plane (see Eq.l24li yields: 

C(v,„-5,„),/ = 5„,/-(v„,,-'j -tobsj. (48) 

From Eqs. l25ll36ll4Tl and|43] we see that there are 5 param- 
eters (b, s, k, d, and r) that describe the whole spectral evolu- 
tion, although they have different values at each stage. Starting 
from basic physical principles, we used boundaries for the dif- 
ferent parameters to avoid unphysical results. These boundaries 
are listed in Table |6] 

The negative values for parameter d stand for the possibility 
for an increase in the Doppler factor, D, with radius. Regarding 
r, negative values correspond to collimation, i.e., a decrease of 
the jet radius, and positive values correspond to an expansion 
process. 
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Table 6. Range for spectral parameters allowed to vary in the fits 
to the observed spectrum for all stages. 



h 


s k 


d 


r 


1 to 2" 


2 to 4 1 to 6 


-2 to 2 


-1 to 1 



if r > 0. 



Table 7. Best fit values for spectral evolution modeling of the 
2006 radio flare in CTA 102, parameters b, d, s, k r and /. 



2005.60-2005.95 2005.95-2006.30 2006.30-2006.80 
ClAl C2 A2 



h 
d 
s 


1 ri+O.OS 

0.2±0.02 
2.1 


1 TIT +0.65 
^■■'-^-0.35 

-0.1±0.03 
2.0 


1.7±0.2 

A 7+0.08 
"•^-0.05 

2.4 


k 


2.7±0.14 


not fitted 


4.7±0.4 


r 


0.60+0.03 


0.35±0.02 


0.90±0.07 


t 


0.02+0.01 




not fitted 



" value hits lower boundary. 



5.1.2. Spectral evolution between 2006.0 and 2006.3 

The second peak in the Sm - Vm plane shows a similar behavior 
to a Compton stage, as stated above. Therefore, we applied the 
equations of the Compton stage to the spectral evolution between 
2006.0 and 2006.3. Since the Compton stage can be explained 
within a 4-dimensional parameter space (note that it does not 
depend on parameter k, Eqs. l25]-l36Tl, we selected carefully a 
physically meaningful combination of parameters. 

Boundaries for the slope of the relativistic electron distribu- 
tion, 5, can be derived by using the evolution of the optically thin 
spectral index of the emission (5=1- 2ao) between 2006.0 and 
2006.3. This evolution is shown in Fig.|6] We obtained values of 
■Smin = 1.4 and imax = 2 from the optically thin spectral indeces 
ao = -0.20 and ao = -0.47, respectively. 

5.1.3. Spectral evolution after 2006.3 

The spectral evolution after 2006.3 shows the typical behavior 
of an adiabatic loss stage. Again, we used the evolution of the 
optically thin spectral index, ao, shown in Fig. |6]to derive lim- 
its for the parameter s. The obtained values are s^i^ = 1.6 and 
■Smax = 2.5. This stage could be divided into two substages (be- 
fore and after 2006.45) but the limited time sampling (with only 
three data points in the first part) forced us to perform the spec- 
tral fitting to the whole adiabatic stage. 



On top of the listed limitations, the evolution of the optically 
thin spectral index, ao, can be us ed to provide estimates of the 
parameter 5=1- Iuq. iMarsch er & Gear (1985) used a lower 
limit of 5 = 2 to keep the shock non-radiative. A flatter spec- 
tral slope, e.g., 5 < 2, would increase the amount of high energy 
electrons, which will dominate the energy density and the pres- 
sure in shock. Since these particles suffer radiative cooling, their 
energy losses would affect the dynamics leading to a radiative 
shock. 

As mentioned before, allowing non-conical jet expansion 
{r + \) also affects the boundaries for the spectral parameters 
k and d presented in Table |6] which depend on the value of r 
in the non-conical case. Note that this should be taken into ac- 
count in order to properly interpret the evolution of the physical 
parameters with distance to the core. 

We developed a least-square algorithm to fit the three differ- 
ent radiation loss stages (Compton, synchrotron, or adiabatic, or 
combinations of them) to the observed evolution using one sin- 
gle set of parameters. During the optimization of x'' observed 
data points were used for calculating the constants Cy^^^. and cj,, ,, 
which allows for an improved This technique allowed us to 
test and analyze different possible scenarios (b, s, k, d, and r). 

5.1.1. Spectral evolution before 2006.0 

The spectral evolution until 2006.0 followed approximately 
the standard evolution described by the shock-in-jet model. 
However, we have not found evidence for a plateau phase of 
turnover flux density, S „, neither after the first Compton stage 
nor after the second Compton-like stage. Therefore, we excluded 
the synchrotron stage from our modeling using a direct transition 
from the Compton to the adiabatic one. The result of this mod- 
eling is presented in Table|2]as model ClAl. 



5.2. Final model and error analysis 

Table |7] lists the best fits and errors of the different stages. In 
the first column, we show the values for ClAl. Since we have 
not found evidence for a synchrotron stage, we assumed a direct 
transition between a Compton and an adiabatic stage. 

A study of t he un certainties of the spectral parameters as in 
Lampto n et al.l (11976') is not suitable due to the small number 
of data points and the strong mathematical interdependence of 
the parameters. This yields to mathematically correct, but non- 
physical solutions. For the same reason, we did not perform an 
analysis of the large 5 or 6-dimensional parameter space. To pro- 
vide first-order error estimates we followed this approach: The 
set of parameters derived minimizes the ;\f^-distribution in the 
parameter space, so we investigated the stability of this point in 
the parameter space. This error analysis is based on the variation 
of the for a given parameter sweep within the listed bound- 
aries (Table |6]l, while keeping the others fixed. The final values 
for the uncertainties were obtained by calculating the 68% prob- 
ability values, assuming a normal distribution for the values of 
the parameters. Note that the distribution is not always sym- 
metric around the minimum value, and this leads to lower and 
upper limits for the error estimates. Fig. Q shows the result of 
our spectral modeling. 

5.3. Modeled spectra and light curves 

With the derived parameters and their error estimates we calcu- 
lated the evolution of flaring spectrum and compared it to the 
observed flux density values at any given epoch. For each time 
step during the flaring activity the turnover frequency and the 
turnover flux density can be thus calculated. The optically thin 
spectral index follows from the assigned spectral slope, s. The 
calculation of the optically thick spectral index, a,, can be per- 
formed following the approach of Tiirler et alj (l2000l) . setting 
~ h/ni (see Eqs. [30l and IZTTi. This leads to a set of four 
spectral parameters (S,,, , v,„ , ao , a,) from which the shape of the 
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Fig. 7. Temporal spectral evolution of the 2006 radio flare in CTA 102 left: turnover frequency; right: turnover flux density. The 
lower panels show the residuum for the fits x = (-^obs - -"Cmodei) / Ax. The dashed black lines correspond to the time labels in Figure 
Eland indicate the extrema in the evolution. 



flaring spectrum can be derived by using EqJT] By incorporat- 
ing the quiescent spectrum, the total spectrum can be computed. 
From the uncertainties of the parameters (b, s, k, d, and r) and 
the quiescent spectrum, the lower and upper boundaries for the 
modeled spectra can be obtained. Fig.|8]shows the modeled spec- 
trum for the 2006.2 observations as an example. The calculation 
of the optically thick spectral index, using a, ~ f^ln-i led to val- 
ues which reproduce well high frequency part of the spectrum 
but not the low frequency one (see Fig|8]l. The discrepancy in 
the optically thick part of the spectrum was probably caused by 
the quiescent contribution, which is known to vary slightly over 
time (see variations in the flux density for t < 1995 Figure|2]l. 



Once the spectra are calculated, it is straightforward to ob- 
tain the modeled light curves at a given frequency. The compar- 
ison between the modeled and the observed 37 GHz light curve 
is presented in Fig. |9] Although similar plots can be obtained 
for other frequencies, it is meaningful to select a frequency at 
which the observed emission is mainly generated by the interac- 
tion of the traveling shock with the underlying flow. Moreover, 
the dense sampling of the 37 GHz observations provided a more 
complete picture of the evolution of the flaring event. Taking this 
consideration into account, we picked the 37 GHz light curve as 
the best possible comparison. The observed data points fall well 
within the range of modeled light curve. At the beginning of 
the flare, the 37 GHz flux density was still in its quiet state and 
therefore the main uncertainty for this time is due to the uncer- 
tainties in the quiescent spectrum. The error band presented here 
corresponds only to the flaring state, which was the dominant 
contribution to the 37 GHz flux density after 2005.60, leading to 
broader error bands from this epoch and on. 



• data 
^— total spectrum 
— - flaring spectrum 

quiescent spectrum 
error spectrum 




Flaring Spectrum: 
t:2006.20 [yr] 

5„:7.76r„"3? [Jy] 
v^:57.nZ^:^^ [GHz 



-0.50 



V [GHz 



Fig. 8. Modeled spectrum for the 2006.20 observations. The 
dashed-dotted line corresponds to the quiescent spectrum, the 
dashed line to the flaring spectrum, and the solid black line to 
the total spectrum. The gray lines indicate the uncertainties in 
the calculation of the spectrum (for more details see Sect. 15.31 
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37 GHz (model) 
37GHz (data) 




0.0' ^ 

2005.60 



2006.20 

t [yr] 



Fig. 9. Modeled and observed 37 GHz light curve. The black 
solid line corresponds to the calculated model and the dashed 
lines indicate the uncertainties in the calculation of the modeled 
light curve (for more details see text). 



5.4. Rejected solutions 

As already mentioned, the solutions for the spectral fitting were 
embedded in a 4- or 5-dimensional parameter space. Some of 
those can be highly degenerate given the small number of data 
points. When fitting the spectral evolution, a deep study of the 
parameter space was performed, including all possible combina- 
tions of the parameters, (i.e., variation of single, pairs and triplets 
of parameters) and simultaneously fitted stages, (one-, two- and 
three-stage fits). Within this study we found families of solu- 
tions with values similar to the presented final set of parame- 
ters, but with unphysical values. After a sanity check, those fits 
were discarded. The meaningful set of solutions is presented in 
Table |2] and discussed in Sect.|6l 



5.4.1. The geometrical model 

IStevens et al.l (Il996 h found a similar double hump in the v„, - S ,„ 
plane for 3C 345. They assumed that this behavior could be due 
to changes of the viewing angle, expressed in a variation of the 
Doppler factor, D, along the jet. Therefore we used as an al- 
ternative to the previous modeling of the second hump in the 
Vm-Sm plane, a purely geometrical approach. Such an approach 
could easily explain the variation in the observed turnover flux 
density (S'„, oc D^^""S,„^ while the observed turnover frequency 
kept a nearly constant value (v',, oc Dv„,). It was assumed that the 
deviation from the standard shock-in-jet model around 2005.95 
was caused by changes in the evolution of the Doppler factor 
cos §) ' during the final adiabatic loss stage. The 
remaining parameters were taken from the best fit to the evolu- 
tion until 2006.0 (first column in Table|7] model CI Al) and only 
the parameter d was allowed to vary. 

The result of our calculations for the different models is pre- 
sented in Fig. [To] and Table [8] The figure shows the fits sepa- 



rating the diff'erent stages observed after 2006.0 (Figs.|4]and|5). 
We recall that these stages are the increase of the flux density 
between 2006.0 and 2006.3 and the decrease in flux density and 
frequency from 2006.3 until 2006.8. 

Table 8. Results of the geometrical model for parameter d. 



time 


type 


d 


2006.0- 


-2006.8 


3-stage 


-0.04 


2006.0- 


-2006.5 


2-stage 


-0.30 


2006.5- 


-2006.8 


2-stage 


-1.24 


2006.0- 


-2006.3 


1 -stage 


-041 


2006.3- 


-2006.5 


1 -stage 


-1.33 


2006.5- 


-2006.8 


1 -stage 


-1.27 



It is obvious that a simple change in the evolution of the 
Doppler factor, expressed by the exponent d, {d oc /? ''^''j, can 
not explain the observed temporal evolution of the turnover fre- 
quency and turnover flux density. 

5.4.2. Applying the modified compton stage 

As mentioned in Sect.|3]lBTomsson & Aslakse"nl ( |2000|) reviewed 
Marscher & Gear (1985) assumption for the Compton stage and 
modified the equations for the evolution of the turnover fre- 
quency, v„,, and the turnover flux density, We applied this 
model to the first two stages of the 2006 flare in CTA 102 (Model 
ClAl), using the expression as in Sect. l3.2l for the time evolu- 
tion (Eq. l39l l and d - Q. The result of the modeling is shown in 
Fig.E] FigHniand Table|9] 
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Fig. 11. The 2006 radio flare in the turnover frequency - turnover 
flux density plane. The time labels indicate the time evolution 
and the temporal position of local and global extrema in the 
spectral evolution. The black line corresponds to the final model 
presented in Tabl^ (based on Marscher & Gear 1985) and the 
red one to Table|9](based on Biomsson & Aslaksenil2000ft . 
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Fig. 10. The temporal evolution of the 20 06 radio flare in CTA 1 02 m odeled with varying Doppler factor (black solid line) for 
/ > 2006.0 and with the modifications of Biornsso n & AslaksenI ('2000) for the first two stages (red solid line). Both models fail 
to describe the 2006 flare (compare to Fig. |7}.; left: turnover frequency; right: turnover flux density. The lower panels show the 
residuum for the fits x - (xohs - J^modei) /Ajc. The dashed black lines correspond to the time labels in Figure |4] and indicate the 
extrema in the evolution. 



Table 9. Spectra l mrameters for ClAl using 

iBiornsson & AslaksenI (l2000l) 



b 


s k 


d r t 


1.0 


1.6 2.9 


0.6 0.02 



The result of th e fitting shows that the approach of 
iBjornsson & Aslaksenl (2000) cannot explain the steep rise of the 
turnover flux density with decreasing turnover frequency in the 
first stage of the flare (see solid red line in Fig. [TTl . Following 
the authors, this could be an indication that inverse Compton 
scattering is not th e dominant energy loss m echanism during the 
rising stage (Bior nsson & AslaksenI l2000b . This aspect will be 
discussed in Sect.|6] 

6. Discussion 

Our analysis of the strong radio flare observed in CTA 102 
around 2006, shows that its behavior until epoch 2006.0, i.e., 
the first hump in the - Sm plane, can be well modeled by 
a Compt on and adiabatic stage within the standard shock-in-jet 
model (Mar scher & Geail[T985l) (see Sect. 15.1. It . We found no 
evidence of a synchrotron stage between the Compton and the 
adiabatic stages during the first part of the evolution. This be- 
havior was similar to the one foun d in another prominent blazar 
3C 345 dLobanov & ZensusI 19991) . where a bright jet component 
also appeared to proceed from the Compton stage to the adia- 
batic stage, with the synchrotron stage ruled out by combination 
of the observed spectral and kinematic evolution of this compo- 
nent. The flare in CTA 102 shows a second hump in the v,„ - S,„ 



plane (see FigfTTI) after 2006.0, that cannot be explained within 
this model. The evolution between 2006.0 and 2006.3 seems to 
follow the predicted evolution of a Compton stage, i.e., increas- 
ing turnover flux density and decreasing the turnover frequency. 
Based on this apparent behavior we applied the equations of 
the Compton stage to this second hump (see Sect. 15.1.2b . After 
epoch 2006.3, the evolution shows a decrease in both turnover 
values, which can be interpreted as an adiabatic stage. Thus, we 
used the appropriate equations for this stage to fit the data points 
(see Sect. 15.1.31 ). This approach led to a set of exponents which 
describe the evolution of the magnetic field, B, the Doppler fac- 
tor, D, the spectral slope of the electron distribution, s, the nor- 
malization coefficient of the relativistic electron distribution, K, 
and the jet radius, R. The results of our modeling are summarized 
in Fig. |7] Our hypothesis is that this behavior can be interpreted 
in terms of the interaction between a traveling shock and a stand- 
ing shock wave in an over-pressured jet. A pure geometrical 
model failed to explain the observed behavior (see Sect. 15.4.1b . 
We also tried to fit the evolution of the flare using the modifica - 
tion of the shock-in-jet model bv lBjornsson & AslaksenI (l2000l) . 
but this also failed. However, the failure of this model, which 
takes into account multiple Compton scattering, could have fur- 
ther implications in our understanding of flaring events in AGN 
jets, as discussed at the end of this section. 

The formation of a standing shock can be described in the 
following way: The unbalance between the jet pressure and the 
pressure of the ambient medium at the jet nozzle leads to an 
opening of the jet. Due to the conservation laws of hydrodynam- 
ics, this opening results in a decrease of the density, the pressure 
and the magnetic field intensity in the jet. The finite speed of 
the sound waves in the jet is responsible for an over-expansion 
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followed by a re-collimation of the jet that gives rise to the for- 
mation of a shock. During this collimation process the jet radius 
decreases and the shock leads to an increase of the pressure, den- 
sity, and magnetic field intensity. Again, the finite speed of the 
sound waves is responsible for a over-collimation of the jet. This 
interplay between over-expansion and over-collimation leads to 
the picture of a pinching flow, i.e., a continuous change of the 
width along the jet axis, in contrast to conical jets. The intrin- 
sic physical parameters (pressure, density, and magnetic field) 
alon g a pinch ing jet show a seque nce of local maxima and min- 
ima dDaiiL&Marscher 1988; FaU"g [T99l . 

Summarizing, an over-pressured jet can be divided into three 
regions, i) the expansion region, i.e., continuous increase of the 
jet radius, ii) the collimation region, i. e., decrease of the jet 
radius and formation of the collimation shock, and iii) the re- 
expansion region (see Fig.fTSli. In such a scenario, enhancements 
of emission can be produced by the interaction between traveling 
and standing (re-collimation) shocks (.Gomez et al.til997,) . In the 
following, the evolution of the travelling shock in these regions 
is described and our results are put in context. 



Pa PO > Pa 




Fig. 12. Sketch of an over- pressured jet with indicat ed character- 
istic regions (adopted from lDalv & Marscherlll988l) . 



6.1. The expansion region 

A relativistic shock propagating through this region accelerates 
particles at the shock front. These particles travel behind the 
shock and suffer different energy loss mechanisms, depending on 
their energy. The resulting evolution of the turnover frequency 
and turnover fl ux density is explained by the shock-in-jet model 
(iMarscher & G ear 1985) under certain assumptions. 

The parameters derived for the time between 2005.6 and 
2006.0 can be associated with this region, including a Compton 
and an adiabatic stage. The expansion of the jet is parametrized 
by r = 0.60, which differs from the value expected for a conical 
jet r = 1. The n on-conical behav ior could be due to accelera- 
tion of the flow (lMarscheiill980l) . This acceleration should be 
expressed hy -r ■ d > {D <x L ''''). However, from our results, 
this value is negative -r ■ d - -0.12, but very small, i.e., com- 
patible with no changes in the Doppler factor Thus, we cannot 
confirm this point. 

For the evolution of the magnetic field with distance we de- 
rive a value of b - 1 .0, indicating that the magnetic field could 
be basically toroidal in this region. The injected spectral slope 
for the relativistic electron distribution s = 2.1 leads to an opti- 
cally thin spectral index oq - -0.55. A decrease in the density 
can be deduced from -r-k = -1.6, which corresponds to the evo- 
lution of the normalization coefficient of the relativistic electron 
distribution, K. 

The parameter foff corresponds to the time difference be- 
tween the onset of the Compton stage and the first detection of 



the ffare. From the value of foff = 0.02 yr together with inde- 
pendently obtained values for the viewing angle, § - 2.6°, and 
the apparent speed of the VLBI compone nt ejected by the 2006 
flare, jS^pp = 17 c dJorstad et all 120051: iFromm et allfeoiO ). we 
calculated the displacement between the onset and the detection 
of the flare to be Ar = 3.5 pc. 



6.2. Tine coiiimation region 

After the recollimation region, at the position of the hypothet- 
ical standing shock, the local increase in density, pressure and 
magnetic field should generate an increase in the emission. The 
interaction between a trav elling and a standing shock would fur- 
ther enhance the emission dGomez et al.ll997 ). Furthermore, the 
standing shock would be dragged downstream by the traveling 
shock and re-established after a certain ti me at its initial position 
(iGomez et al ] ll997l:lMimica et alJl2009h . We compare here our 
results with this scenario. 

Since the evolution of the turnover frequency and turnover 
flux density between 2006.0 and 2006.3 showed Compton-stage- 
like behavior, i.e., decreasing turnover frequency and increasing 
flux density, we used the equations of the Compton stage to de- 
rive the possible evolution of the physical parameters (model 
C2). In this region a slower rate of jet expansion is found 
(r = 0.35). During this stage, the Doppler factor seems to be 
constant with distance, -r ■ d - 0.035. In the context of the 
hypothetical shock-shock interaction, acceleration of the flow 
close to the axis is expected down to the discontinuity of the 
stading shock, where su dden deceleration would occur (see, e.g., 
IPerucho & Martd2007h . Thus, it is difficult to assess whether the 
Doppler factor should increase or decrease in the whole region. 

The magnetic field intensity decreases with an exponent 
b - 1.35, implying that the geometry of the magnetic field 
has changed, with contributions of non-toroidal components, but 
showing no hints of magnetic field enhancement. The parameter 
s, giving the spectral slope of the relativistic electron distribu- 
tion changes to i = 2, which gives an optically thin spectral 
index ao = -0.5. 

The set of parameters derived for this time interval do not 
reflect the expected physical conditions of a traveling -standing 
shock interaction, other than a slight flattening of the spectral 
slope (from possible refreshment of particles). Nevertheless, the 
shock-shock scenario could hardly be reproduced by a one- 
dimensional model. Numerical simulations should be performed 
in order to study this hypothesis in detail. Another possibility is 
that the reason for the second peak is attached to the injection of 
a second shock from the basis of the jet. This is not observed, 
though. 



6.3. Tine re-expansion region 

After the re-collimation process, the jet re-expands, i.e., the 
jet radius increases again. In principle, the position of the re- 
collimation shock can be regarded as a "new" nozzle from which 
the fluid emerges. Therefore, when the shock front reaches this 
region, the expected evolution is, again, that predicted by the 
shock-in-jet model. 

The evolution between 2006.3 and 2006.8 is identified 
within our hypothesis with the re-expansion region. Thus, the 
equations for an adiabatic loss stage were applied to the evolu- 
tion turnover frequency and turnover flux density. 

The opening of the jet is clearly apparent at this stage r = 
0.90. This opening should produce a decrease in density, which 
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translated into smaller values of the parameter -r - A;, cc L j. 
From the fits, we derive -r - k - -4.2, confirming a decay in the 
density. The magnetic field falls with b - 1 .7, which shows again 
that the geometry of the field changes from a purely toroidal 
to a mixed structure with the distance. The values for -r ■ d - 
0.18 reveal an acceleration of the flow, which can naturally arise 
during the expansion of the jet. The spectral slope of i = 2.4 
translates into an optically thin spectral index of ao - -0.7. 

6.4. Spectral slopes and optically thin spectral indices 

The variation of the fitted optically thin spectral index, orof, 
(see upper panel in Fig. |6]l is an indication of the aging of 
the relativistic electron distribution due to the diff'erent energy 
loss mechanisms during the evolution of the flare. The model 
presented by Marscher & Gear (1985 ) did not include such an 
aging of the relativistic electron distribution. Despite this re- 
striction, we could model the qualitative behavior of the evo- 
lution via the different values obtained for parameter s in the 
different fitted stages: For the expansion region we obtained 
ao.m - -0.55, a slightly flatter spectral index for the coUimation 
region ao,m = -0.50 and a steeper value for the re-expansion 
ao,m = -0.7. 

6.5. The influence of the interpolation and the quiescent 
spectrum 

We applied several interpolation steps together with different 
quiescent spectral parameters to the analysis of the light curves 
in order to test their influence on the study of the evolution 
of the peak parameters in the turnover frequency -turnover flux 
density (v,„ - 5,,,) plane. All the tests showed a second hump 
in the v,„ - 5,„ plane. From the tests, the positional shifts (in 
time of appearance, turnover frequency and turnover flux den- 
sity) of the second hump were calculated: The hump appeared 
att- 2006.30 + 0.05 yr at a turnover frequency v,„ = 56;^3 GHz 
and turnover flux density 5,„ = 8.5^g g Jy. 

These tests proved that the increase of the turnover frequency 
and turnover flux density around 2006.3 is not an artifact gener- 
ated by the interpolation and/or the choice of the quiescent spec- 
trum. 

6.6. Adiabatic versus non-adiabatic shock 

iBiornsson & AslaksenI (|2000|) suggested that the steep slopes in 
the rising region of the Vm-Sm- plane, which is generally iden- 
tified as the Compton stage in the shock-in-jet model, could 
not be due to Compton radiation, but to other processes such 
as isotropization of the electron distribution or even to the non- 
adiabatic nature of the shock. 

The slope we obtained for this time-interval in the v,„-5,„ 
plane is too s teep to be reproduced by the model presented in 
IBiornsson & Aslaksen (2000) (see Sect. 15.4.^ . Thus, the non- 
adiabatic nature of the first stages of evolution of shocks in ex- 
tragalactic jets remains an open issue on the basis of our limited 
data set. 

7. Summary and conclusions 

In this work we present the analysis of the 2006 flare in CTA 102 
in the cm - mm regime. The obtained evolution could be well 
described, up to a c ertain point (t < 2006.0 ), with the standard 
shock-in-jet model dMarscher & Gear|ll985l) . In order to model 



the further evolution of the flare we proposed a second Compton- 
like stage and a second adiabatic stage. We derived the evolution 
of the physical parameters of the jet and the flare and we per- 
formed a parameter space analysis to obtain the uncertainties of 
the values obtained. From the modeled parameters, together with 
their uncertainties, the theoretical light curves and spectra were 
computed. The result was shown to be in fair agreement with the 
observations. 

However, the shock-in-jet model is not able to reproduce the 
second peak of the double-hump structure found in the evolu- 
tion of the peak flux - peak frequency plane in a consistent way 
within our hypothesis of a shock-shock interaction. This is surely 
due to the complexity of the situation and does not rule out our 
hypothesis. Numerical simulations will be performed to test it. 
The goodness of the fit in this region may be due basically to the 
lack of observational data as compared to the number of param- 
eters implied in the modeling. 

Moreover, our results show that some of the basic assump- 
tions of the present models should be reviewed. In particular, 
we found that the first stage of the studied evolution could 
be incompatible with an ad iabatic behavior if the model of 
IBiornsson & AslaksenI (|2000|) is correct. 

We plan to improve our results by performing a similar anal- 
ysis to multifrequency VLBI observations. Future work on the 
spectral evolution of components in jets, especially the interac- 
tion between traveling shocks and re-collimation shocks, should 
include numerical simulations to help to understand this non- 
linear phenomenon, and improve the strong limitations of the 
assumptions required by the shock-in-jet and other analytical 
models. The analysis of multifrequency VLBI observations dur- 
ing the 2006 radio flare in CTA 102 will be presented in paper 
II. 
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